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be ABSTRACT 

■ We present results of an observation with Suzaku of the dipping, periodic bursting low mass X-ray binary XB 1323-619. Using 
the energy band 0.8 - 70 keV, we show that the source spectrum is well-described as the emission of an extended accretion disk 
^ . corona, plus a small contribution of blackbody emission from the neutron star. The dip spectrum is well-fitted by the progressive 
QQ covering model in which the extended ADC is progressively overlapped by the absorbing bulge of low ionization state in the 
^-H outer accretion disk and that dipping is basically due to photoelectric absorption in the bulge. An energy-independent decrease 
' of flux at high energies (20 - 70 keV) is shown to be consistent with the level of Thomson scattering expected in the bulge. 
' An absorption feature at 6.67 keV (Fe XXV) is detected in the non-dip spectrum and other possible weak features. In dipping, 
1^-^ , absorption lines of medium and highly ionized states are seen suggestive of absorption in the ADC but there is no evidence 
(^*^ ■ that the lines are stronger than in non-dip. We show that the luminosity of the source has changed substantially since the 
' Exosat observation of 1985, increasing in luminosity between 1985 and 2003, then in 2003 - 2007 falling to the initial low value. 
. X-ray bursting has again become periodic, which it ceased to do in its highest luminosity state, and we find that the X-ray 
^ ' bursts exhibit both the fast decay and later slow decay characteristic of the rp burning process. We present arguments against 
. ^ , the recent proposal that the decrease of continuum flux in the dipping LMXB in general can be explained as absorption in an 
1 ionized absorber rather than in the bulge in the outer disk generally accepted to be the site of absorption. 

. . . I Key words. Accretion: accretion discs - acceleration of particles - binaries: close - line: formation ~ stars: neutron - X-rays: 
binaries - X-rays: individual (XB 1323-619) 



1. Introduction 

XB 1323-619 is a member of the dipping class of Low Mass 
X-ray Binaries (LMXB), exhibiting reductions in X-ray in- 
tensity at the orbital period due to absorption in the bulge 
of the outer accretion disk. (White & Swank 1982; Walter 
et al. 1982). It is remarkable as one of the two "clocked" 
X-ray bursts sources in which bursting remains close to 
periodic over extended periods of time, the other source 
being 4U 1826-24 (Galloway et al. 2004). Other sources 
such as 4U 1705-44 (Langmeier et al. 1987) have shown 
bursting close to periodic, but this is not sustained over 
long periods of time (Cornelisse et al. 2003). The dip- 
ping sources display relatively complex spectral evolution 
in dipping, which provides strong constraints on emission 
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models for LMXB in general since the model must be able 
to fit the non-dip state and a sequence of spectra formed 
by selection in intensity bands at different depths of dip- 
ping. Work over a period of 15 years provided substantial 
evidence for a model in which the X-ray emission con- 
sists primarily of Comptonized emission from the accre- 
tion disk corona (ADC) plus simple blackbody emission 
from the neutron star (Church et al. 1997, 1998a, b, 2005; 
Balucihska-Church et al. 1999, 2000a; Smale et al. 2001; 
Barnard et al. 2001). This work establishes that in the dip 
sources in general, the dominant Comptonized emission is 
removed slowly showing that the emission is extended, 
while the point-like blackbody emission which we asso- 
ciate with the neutron star is rapidly removed, and causes 
rapid variability in dipping. Measurement of dip ingress 
times allows the size of the extended emitter to be mea- 
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sured , and this has shown the ADC to reach over a large 
part of the inner accretion disk, have a radial extent R be- 
tween 20 000 and 700 000 km, and is thin, with height H 
having H/R « 1 (Church & Balucinska-Church 2004). 
Strong independent support for the extended ADC comes 
from the Chandra grating results of Schulz et al. (2009) in 
which emission lines of highly ionizied species were found 
to originate at radial distances between 20 000 and 110000 
km in an extended ADC. Various authors have addressed 
theoretically the question of the formation, size and prop- 
erties of an ADC, such as Jimencz-Garatc et al. (2002, 
and references therein) by considering the effects of illu- 
mination by the central source. 

In the present work we investigate spectral evolution in 
dipping utilising the broadband capability of Suzaku, the 
long-term luminosity evolution of the source and the asso- 
ciated strong changes in burst properties. In general, the 
X-ray burst sources do not have regular bursting. If it is as- 
sumed that the mass accretion rate does not vary greatly 
on timescales of a few hours, regular bursting would be 
expected, since mass accumulates at a steady rate on the 
surface of the neutron star and conditions for unstable 
nuclear burning arc achieved. However, the majority of 
sources have quite irregular bursting and the study of the 
two "clocked bursters" is clearly of interest in explaining 
this. The source was observed by us using XMM-Newton in 
2003, January, and analysis of this observation and com- 
parison with observations previously made with Exosat, 
ASCA, and three observations with RXTE revealed rela- 
tively large systematic changes in burst properties and in 
the luminosity of XB 1323-619 over the period 1985 - 2003 
(Church et al. 2005). The most striking feature was the 
changing burst rate, the time between bursts At decreas- 
ing linearly with time from 322 minutes in 1985 to 58 min 
in 2003. During this period the X-ray luminosity in the 
band 1-10 keV increased substantially from 1.3 x 10'^^ 
to 5.4 X 10^^ erg s~^. The observed trend if continued 
would mean that the burst rate would continue to increase, 
implying that At would become zero on Jan 11, 2008. 
However, analysis of a further observation with RXTE in 
2003, September showed a sudden change in the source, 
the luminosity decreasing and the time between bursts in- 
creasing to a mean of 162 min (Balucinska-Church et al. 
2008). 

XB 1323-619 was discovered using Uhuru and Ariel V 
(Forman et al. 1978; Warwick et al. 1981) and X-ray 
bursting and dipping found using Exosat by van der 
Klis et al. (1985). An orbital period of 2.938± 0.020 hr 
was determined by Balucinska-Church et al. (1999) in 
good agreement with the value of 2.932± 0.005 hr from 
Exosat (Parmar et al. 1989). The source was observed with 
BeppoSAX in 1997 from which a detailed study of the 
broadband 1 - 200 kcV spectrum was made (Balucinska- 
Church et al. 1999) showing the Comptonized emission 
to have a relatively high cut-off energy of 44 keV. The 
XMM observation (Church et al. 2005) revealed various 
line and possible edge features, including Fe absorption 
lines at 6.70 and 6.98 keV equivalent to Fe XXV and 



Fe XXVI K„ and an emission line at 6.82 keV. In addi- 
tion absorption features were seen at 1.46 keV and 4.05 
keV, possibly Mg XII and Ca XX. The origin of lines was 
discussed and it was concluded that the evidence favoured 
production in the ADC. 

In analysis of the XMM data, Boirin et al. (2005) and 
Diaz-Trigo (2006) suggested that the dipping in all LMXB 
dipping sources could be explained in terms of an ionized 
absorber rather than in the low ionization state bulge in 
the outer disk generally accepted as the site of absorption. 
This model woiild require dipping to be unconnected with 
the bulge in the outer disk. In the present observation 
we show that dipping is very well explained as absorption 
in the bulge plus a small degree of energy-independent 
electron scattering consistent with the measured column 
density of the cool absorber. In Sect. 4 we discuss general 
objections to the ionized absorber model. 

2. Observations and analysis 

We observed XB 1323-619 with Suzaku in 2007, Jan 9 - 10 

for 34.1 hours, this covering 11.7 orbital cycles. Data from 
both the XIS and HXD instruments were used. Analysis 
was carried out using the Suzaku-speci&c Ftools, part of 
FTOOLS 6.3.1. At IS AS, the XIS raw telemetry files were 
converted to Fits format using the facility mkIstfits, and 
then the raw data files were converted to filtered and cal- 
ibrated events files using mk2ndfits. This second stage 
was carried using the improvements in calibration data of 
August, 2007 (version 2 processing). The XIS was oper- 
ated in the normal mode, using the one-quarter window 
option viewing 1/4 of the CCD and with the 3x3 and 5x5 
pixel options for pulse height analysis and with 2 s time 
resolution. Spaced-row charge injection (SCI) was used to 
reduce the effects of radiation damage. Data for each of 
the three detectors XISO, XISl and XIS3 detectors were 
available, XIS2 having failed prior to this observation. The 
filtering consisted of the removal of hot and flickering pix- 
els and of selection on event "grade" from the pulse height 
distribution (keeping grades 0, 2, 3, 4 and 6) to discrimi- 
nate between X-ray and charged particle events, and also 
standard screening based on housekeeping data as follows. 

This standard screening excluded periods of South 
Atlantic Anomaly passage plus an elapsed time of 436 
seconds after the time of leaving SAA (T_SAA > 436), 
to select data for elevation above the Earth's limb more 
than 5 degrees (ELV > 5), for elevation above the sunlit 
Earth more than 20 degrees (DYE_ELV > 20), and for ge- 
omagnetic rigidity (COR > 6) . Lightcurvcs were extracted 
from each XIS detector from the filtered events file using 
Xselect. The lightcurve obtained after screening revealed 
only a single burst and 9 dips, whereas a lightcurve made 
with relaxed screening criteria (below) showed 5 bursts, 
the last of which was double, and eleven dips, so clearly 
standard screening removed substantial amounts of data 
of interest and we investigated the effects of relaxing the 
selection criteria. 
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To do this, it was necessary to use the unfiltered events 
file. Firstly, hot and flickering pixels were removed using 
the FTOOL CLEANSIS and grade selection was made within 
Xselect. Inspection of the housekeeping data showed that 
during the observation, elevation of the source above the 
Earth's limb was always above 7.5°, i.e. no Earth occulta- 
tion took place and so there was no need to screen for this. 
However, wc found that when we reduced the daytime ele- 
vation selection criterion to DYE_ELV > 10 ° , more bursts 
were revealed and we eventually made no selection on this 
parameter. Wc then checked the effect of this on spectral 
data, since X-rays from the bright Earth have a very soft 
spectrum dominated by O and N lines below 1 keV. Non- 
burst spectra with and without screening for DYE_ELV 
were extracted, but no differences between them could 
be seen, particularly below 1 keV. As the spectrum of 
XB 1323-619 in the XIS does not extend below -0.8 keV 
because of its low Galactic latitude source and associated 
large column density, little effect would be expected. 

The geomagnetic rigidity (COR) was also examined. 
With the standard selection: COR > 6, a small fraction of 
data was removed including the first burst at time ~6,000 
s (sec Fig.l) when the COR value was —5. Therefore to 
include all the bursts we relaxed the selection criterion 
to COR > 5. This may slightly increase the background, 
however as XB 1323-619 is a bright source in XIS, the over- 
all effect will be negligible. It is clear that our approach 
of relaxing the selection criteria will not affect the deter- 
mination from the lightcurvc of the burst repetition rate 
and decay times, and our testing showed that the effect 
on spectra was negligible. 

With the final selections above, the lightcurves from 
each XIS detector were examined for traces of the 
nearby pulsar ISAX J1324.4-6200 (Angelini et al. 1998; 
Balucihska-Church et al. 1999) located 17'from XB 1323- 
619 which would have contaminate lightcurves and spec- 
tra. The observation had been designed with a pointing 
direction to minimize possible effects. No trace of pul- 
sations could be seen. A lightcurve of XB 1323-619 was 
then extracted using all 3 XIS detectors which were each 
operated in the mode in which one quarter of the full 
detector was used consisting of a strip of 1024 x 256 pix- 
els. The lightcurve was extracted in the full energy band 
0.2 - 12 keV by specifying an approximately rectangular 
source region with the source at the centre of size 4.4'x 
6.7'; background was extracted using two square regions 
in the image, each side of the source at a distance of 7' 
and 4.4'square, and adding these. Background subtraction 
was carried out using lcmath allowing for the different 
areas of the somcc and background regions. Photon pileup 
is negligible. The background-subtracted lightcurve of the 
3 XIS detectors is shown in Fig. 1. 

HXD data were obtained in the form of a cleaned 
events file in which standard screening was carried out 
for SAA passage, elevation etc. We analysed the PIN data 
and found that the source is detected up to ^^50 keV. 
Because of this and the weakness of the source, we did not 
use the high energy GSO data in the band 30 - 600 keV. 



Non X-ray background in HXD/PIN is evaluated using 
a background model provided for observers for each ob- 
servation. The cosmic X-ray background which is about 
5% of the total is included as a model term in spectral 
fitting using a power law form corresponding to previ- 
ous determinations of the spectrum. Data were selected 
via a good time interval file (GTI) ensuring that both 
source and background data were available. Lightcurves 
were background-subtracted using the instrumental back- 
ground, and deadtime corrected using hxfdtcor. DYE_ 
filtering is not carried out and so does not affect the num- 
ber of bursts seen; however, COR > 6 is chosen, so the first 
burst seen in XIS is removed, but the remaining bursts are 
seen. 

3. Results 

3.1. The burst rate 

The lightcurve of Fig. 1 shows 5 bursts, the last of which is 
double. We first obtain the mean time between bursts At 
and thus the burst rate {At~^). This was done by exam- 
ining the XISO lightcurve with its full 2 s resolution, i.e. 
without 64 s binning. The times of the peaks of the bursts 
and separations between the bursts are listed in Table 1. 
The gap between the 1st and 2nd bursts is consistent with 
being twice the approximate separation of bursts, and it is 
almost certain that one burst between these fell in a data 
gap and so wc use half the gap between these bursts. We 
exclude the double burst from the analysis and using the 
gaps between the other bursts we find that At = 20254 ± 
106 s (with 1(7 errors), i.e. At = 5.63 ±0.03 hr. 

Figure 2 (left panel) shows the evolution of burst re- 
currence time At with time, i.e. MJD spanning the period 
1985 to the Suzaku observation. Analysis producing At 
is taken from Balucihska-Church et al. (2008) including 
values derived for the observations with Exosat, ASCA, 
RXTE, XMM-Newton and added to it is the point for the 
present Suzaku observation. It can be seen that the burst 
rate shows a systematic linear decrease from the Exosat 
observation to the XMM observation (2003, Jan 29) as 
shown by the best-fit solid line. The observations of XTE 
2003 (2003, Sep 25) and the present observation do not 
follow the trend as marked with the dotted line. It ap- 
pears that a marked change took place within the source 
in 2003. 

From spectral fitting below using the best fit to the 
broadband spectrum of the non-burst, non-dip data we 
obtain a 1 - 10 keV luminosity of 1.47 x 10^^ erg s~^, and 
we next compare the Suzaku data with our previous work 
on XB 1323-619 in which we obtained the dependence of 
burst rate on source luminosity from the Exosat observa- 
tion of 1985 to the RXTE observation of 2003, a span of 
18 years (Baluciiiska-Church et al. 2008). Burst rate as 
a function of luminosity is shown in Fig. 2 (right panel) 
including our analysis of the observations with Exosat, 
ASCA, RXTE, XMM-Newton and the present Suzaku ob- 
servation. The solid line represents a linear fit to the data 
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Fig. 1. Background-subtracted lightcurve of the January 2007 observation of XB 1323-619in the energy band 0.2 - 12 keV with 
64 s binning from the XIS detectors XISO, XISl and XIS3 
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Fig. 2. Left panel: recurrence time At of X-ray bursts showing the systematic decrease of At from the Exosat observation to the 
XMM observation followed by the departures from that pattern of the 2003 RXTE and the present observations; right panel: 
the dependence of burst rate 1/At as a function of luminosity. There is an approximately linear relation from the Exosat to 
the XMM observation, although the 2003 RXTE and Suzaku points have reverted to the positions consistent with the much 
decreased luminosities 



up to the 2003 XTE observations. Firstly, it can be seen 
that the 2003 XTE point follows the linear relation be- 
tween the luminosity and the burst rate. Moreover, the 
point for Suzaku is close to the Exosat point, indicating 
not only that the source has jumped back to its luminos- 
ity of 1985, but also that the burst rate is consistent with 
that luminosity. 

3.2. The burst properties 

We determined the times of the burst peaks, the time be- 
tween bursts and examined the decays of the bursts and 



for this, the 64 s binning of Fig. 1 was not adequate and so 
we produced a lightcurve from XISO with 2 s binning, this 
having more than sufficient count statistics. The decay 
profiles of the bursts were examined, and it was immedi- 
ately obvious that each burst had a fast decay followed 
by a much slower decay giving each burst a duration of 
more than 200 s. This behaviour was typical of the source 
in previous observations (Balucinska-Church et al. 2008) 
and the slow decay indicates nuclear burning by the the 
rp- process (Wallace & Woosley 1981; Galloway et al. 2004) 
was operating. We fitted each burst separately with a func- 
tion f{t) = A + Ni X exp[-{t - io)/Ti] + x ex'p\-{t - 
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*o)/t2 in the vicinity of each burst, i.e. extending 600 s 
beyond the burst peak. The function contains two time 
constants Ti and t^, to is the burst peak time, A is the 
persistent emission count rate and Ni, N2 are normaliza- 
tions. The results are shown in Tabic 1, from which we 
find that the mean exponential decay times are n = 9 
±1.8 and T2 = 33±7 seconds. 



Table 1. Measured burst properties 



Burst No. 


Peak 
Time (s) 


Separation 
At (s) 


Tl 

(s) 


T2 

(s) 


1 


4148.0 






26.711:1 


2 


44358.0 


40,210.0 


10.4 ±2.8 


97 q+3-9 

^' -^-8.9 


3 


64696.0 


20,338.0 


6 7+2-9 
"■'-2.0 


00 r,+4A 

Zo.Z_2 Y 


4 


85016.0 


20,320.0 


11.1 ± 1.2 


46 1+2* 
^"■'■-12.2 


5a 


107420.0 


22,404.0 


8.61^-2 


35.5^3 Q 


5b 


108054.0 


634.0 


10.8 ± 0.9 


(2) 



times relative to MJD 54109.48056. 
'2) No long decay was present in burst 5b 



The double burst is shown in the Table as bursts 5a 
and 5b; it is interesting that the second weaker component 
5b displayed no evidence for a slow decay. 

3.3. X-ray dipping 

The X-ray dipping was examined by first creating a 
lightcurve from which all portions of all bursts were ex- 
cluded to prevent this interfering with the determina- 
tion of the best-fit orbital period using the dipping. The 
XIS lightcurve which was a sum of background-subtracted 
lightcurves from XISO, XISl, XIS3 in the full band 0.2 - 
12.0 keV was searched for periodicities using the PERIOD 
program which gave a best fit period of 2.926 lo;25 
To reduce the uncertainties, a second method was used in 
which the time of each dip minimum was obtained was 
plotted against cycle number and these data fitted. In the 
present observation, the dips consisted of a first minimum 
followed by a second minimum within the envelope of the 
full dip, so the lowest point in the second part was used, 
there being 8 dips in which the time of this was well- 
defined. This procedure produced a period of 2.928±0.002 
hr. This is in good agreement with the previous determi- 
nations of 2.938±0.020 hr from BeppoSAX (Balucinska- 
Church et al. 1999) and from Exosat of 2.932±0.005 hr 
(Parmar et al. 1989). The data could then be folded on 
this period and in Fig. 3 we show the XIS lightcurves in 
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Fig. 3. Background subtracted XIS lightcurves (the sum of 
XISO, XISl and XIS3) with bursts removed folded on the de- 
rived orbital period of 2.928 hr using 64 phasebins per cycle: 
upper panel: in the band 1.0 - 4.0 keV; middle panel: in the 
band 4.0 - 12.0 keV. The lower panel is the hardness ratio 
formed by dividing the lightcurves in the higher and lower en- 
ergy bands 



two bands: 1.0 - 4.0 keV and 4.0 - 12.0 keV and the hard- 
ness ratio formed by dividing these. The lightcurve was 
folded at the epoch of 54109.577080 MJD, which is the 
time of the first deep dip. Dipping is 55% deep in the lower 
energy band and 36% deep in the higher energy band in 
the folded light curve. However, folding introduces some 
smoothing and in the normal lightcurve (Fig. 1) dipping 
is of the order of 65% deep. The hardness ratio in Fig. 2 
shows that lower energies are removed in dipping leading 
to the observed hardening of the spectrum, which is quite 
typical for the source. 

Similarly, PIN lightcurves extracted in various bands 
were examined and weak dipping with a depth of 10% ± 
6 % can be detected in a folded lightcurve at energies up 
to 20 keV. At higher energies dipping was insignificant. 

3.4. The non-dip non-burst spectrum 

A non-dip, non-burst spectrum was accumulated from the 
whole observation by firstly noting the times of the start 
and end of bursts and removing the bursts by selecting 
on time. To remove dipping intervals and ensure that the 
shoulders of dipping were completely removed a selection 
on phase was made, phases -0.28 to ±0.28 were removed 
where the zero of phase was defined as the deepest point 
of dipping. A spectrum was selected from each of the three 
XIS detectors; the spectrum in each detector has an ex- 
posure time of about 34 ksec giving about 84000 counts 
in each detector (a count rate of 2.5 count s^^). A back- 
ground spectrum was accumulated with the same selec- 
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tions for each of the XIS instruments. A response file was 
generated using the ftool xisrmfgen and a correspond- 
ing auxihary response file also generated for each detector, 
using the latest calibration database of 18th April, 2008. 
The ASCA ftool ADDASCASPEC was then used to add the 
three spectra, the three corresponding background files, 
the three auxiliary response (arf) files and the three re- 
sponse matrix (rmf) files. 

The PIN spectrum was produced using the same GTI 
files as used to select the XIS data. A burst-like event 
was noted in the background subtracted PIN lightcurve 
which occurred ~900 s before burst no. 2 and lasted ~18 s. 
No increase of the count rate was detected in any of the 
three XIS instruments at this time. However, the time 
interval covering this event was removed from the selection 
of the PIN persistent emission spectrum. A background 
spectrum was extracted from the non X-ray background 
events file provided for the observation, using the same 
GTI as above. The cosmic X-ray background was added 
as a term to the spectral model, with photon index F = 
1.29, a high energy cut-off at 40 keV and normalization 
of 8.8 X 10~* ph cm~^ s~^ keV~^ as required for the 
observation with the source at the XIS aimpoint. 

Deadtime correction of the source spectra was per- 
formed using HXDDTCOR. However, this ftool works cor- 
rectly only good time intervals longer than 128 seconds. 
Our PIN spectrum which was selected on orbital phase 
and for removal of bursts had a tendency of having GTIs 
in quite short blocks of time, e.g. 16 s. This problems was 
solved using the standard technique for correcting PIN 
lightcurves, i.e. by making a lightcurve with 128 s bin- 
ning, selecting parts of this with GTIs required and dead- 
time corrected by extracting from the "pseudo event file" 
provided. Systematic errors of 3% were added to all PIN 
spectra; no errors were added to XIS as is normal practice. 

The XIS spectrum was fitted simultaneously with the 
PIN spectrum, parameter values being chained between 
the two. A multiplying factor was included in the model 
to allow for uncertainty in the relative instrument calibra- 
tions and was set to unity for the XIS and free for PIN 
as is usual practice. Channels below 0.8 keV and above 
12 keV in the XIS were ignored where the instrument re- 
sponse is not so well-determined; grouping to a minimum 
of 100 counts per channel was applied. In the PIN chan- 
nels between 12 and 70 keV were used, the instrument 
response at 10 - 12 keV being less certain and above 70 
keV the background was becoming large. 

A number of models were tried. Firstly, a cut-off power 
law gave an acceptable fit with x^/d.o.f. = 1415/1633. 
However, since the 1990s it has generally been accepted 
that LMXB require a non-thermal and a thermal emis- 
sion component as published evidence is very strong. 
White et al. (1988) demonstrated the dominance of the 
Comptonized emission. Curvature in the continuum spec- 
tra at a few keV cannot be explained as a Comptonization 
high energy cut-off and requires the presence of a black- 
body (Church & Balucinska-Church 2001), and as said 
in Sect. 1, the explanation of rapid variability in the 



dip sources in dipping supports the presence of neu- 
tron star blackbody emission. Two-component models 
are clearly needed to fit the most luminous LMXB (e.g. 
Schulz et al. 1989; Hasinger et al. 1990) and fitting two- 
component models became general (e.g. Barret et al. 
2000). Consequently, we next added a simple blackbody to 
the spectral model although in faint sources like XB 1323- 
619 the thermal component is weak. This two-component 
continuum model consisting of an absorbed blackbody and 
cut-off power law ( AG* (cut + bb)) gave a significantly 
better fit for which an F-test showed that it was bet- 
ter with 99.9% probability. Moreover, the normalization 
of the blackbody term was found to be significant at 2.6 
(T, SO from this point we will assume the two-component 
model. However, residuals were present indicating possible 
emission and several absorption features. Because of the 
decreased luminosity of the source compared with that at 
the time of the XMM observation, less line features were 
positively detected. There was some indication of a weak, 
broad emission line at ^6.4 keV; however as discussed be- 
low, we regard detection of this as only tentative. The 
residuals also indicated an absorption feature at ~6.7 keV 
and possible absorption features at ^6.6 and ^6.9 keV as 
shown in Fig. 5. Finally, a trace of a possible absorption 
at ~9 keV can be seen in Fig. 5. 

Before adding lines to the spectral model we re-fitted 
the continuum model excluding the energy range contain- 
ing the lines 6.5 - 8.0 keV to prevent the lines slightly 
distorting the continuum fit. The continuum parameters 
were then frozen and lines added. The width of the emis- 
sion line was fixed (at 0.02 keV) to prevent the line ab- 
sorbing neighbouring continuum. Although the emission 
feature at ~6.4 keV had an equivalent width EW of 17±7 
eV, we do not include it in Table 2 showing the fit results 
as we regard the detection as tentative as the appear- 
ance of the line in the residuals is not totally convincing. 
Initially a single absorption line was added with a fixed at 
0.01 keV as all the absorption features appeared narrow. 
This fitting resulted in a line at 6. 6 (7) ±0.05 keV being de- 
tected with EW = 12±7 eV, and these results are shown 
in Table 2. This line can be identified with Fe XXV and 
so represents absorption by highly ionized material which 
could be located in the ADC. 

However, the residuals indicate there are probably at 
least three weak absorption lines and we attempted to fit 
the spectrum with three lines. In this case, it was not pos- 
sible to determine the energies by free fitting so these were 
fixed at values obtained by inspection of the residuals, i.e. 
6.59, 6.74 keV and 6.90 keV. It was not possible to get 
a sensible fit to the residual feature at 6.9 keV, i.e. the 
normalization was zero. The other two lines at ~6.59 and 
~6.74 keV were found to have EWs of -6.8±6.8 eV and 
-6.4^7 4 eV and correspond to marginal la detections, and 
for this reason we do not show them in Table 2. As said 
above, the energies are uncertain and we should be careful 
not to over-interpret them. The feature at 6.59 keV may 
correspond to a less ionized absorber, that is, to an Fe K 
transition above about Fe XVII. As we see Fe XXV (Table 



Balucinska-Church et al.: Suzaku observations of XB 1323-619 



Suzaku XIS+PDJ XB 1323-619 




Energy (keV) 

Fig. 4. Best fit to non-dip, non-burst spectrum of the three 
XIS detectors fitted simultaneously with the PIN detector; for 
details of the line features detected see Fig. 5 and the text 



2) and also absorption by neutral absorber in dipping, ab- 
sorption by an intermediate ionization state may be quite 
possible. In Fig. 4 wo show the folded data and residuals 
for the full band of the XIS and PIN instruments, while 
in Fig. 5 (upper panel) we show the absorption features 
in the residuals over the band 5-10 keV before lines were 
added to the spectral model. 

In the best-fit solution, the best-fit power law photon 
index F of 1.67^0.11 

is very similar to that we obtained in 
the XMM observation of 1.68±0.08 (Church ct al. 2005) 
and is physically reasonable; the cut-off energy of 85 keV 
is high as expected in a low luminosity source (a value of 
44.ll4;4 was obtained in the broadband observation 
(Baluciriska-Church et al. 1999). In the analysis of the 
XMM observation, absorption at 6.70±0.03 keV was seen 
and 6.98±0.04 keV, and broad emission at about 6.8 keV 
(Church et al. 2005). Barnard et al. (2001) detected an 
emission line at 6.43±0.21 keV in a RXTE observation of 
XB 1323-619. 

3.5. The broadband deep dip spectrum 

It is, of course, of interest to investigate the spectrum of 
deep dipping, to use the broadband capacity of Suzaku to 
reveal the energy dependence of dipping, and test whether 
we can model the spectral evolution. In particular, we ap- 
ply the "progressive covering" model for dipping (Church 
et al. 1997) in which the extended absorbing bulge on the 
outer accretion disk progressively overlaps the extended 
ADC source of Comptonized emission, this model having 
previously provided a very good description of spectral 
evolution in dipping in all dipping LMXB (Church et al. 
1997). We will also test the extent to which dipping can be 
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Suzalsu XIS + PIN XB 1333—610 




5 10 
Energy (IceV) 



Fig. 5. Upper panel: persistent emission spectrum and model 
fit in the band 5-10 kcV before addition of any line fea- 
tures showing several absorption features; lower panel: the cor- 
responding residuals for the deep dip spectrum 



described as photoelectric absorption, and the degree to 
which electron scattering may take place, given the sugges- 
tion by Boirin et al. (2005) based on analysis of the XMM 
observation of XB 1323-619 that spectral evolution in dip- 
ping in all the LMXB dipping sources might be explained 
by photoelectric absorption plus electron scattering in a 
medium of varying ionization parameter. 

Dip data were selected from the deepest part of dipping 
where dipping is ^65% deep, by selecting for intensities 
below 1.2 count in XISO and applying the resulting 
GTI to all instruments. In addition, we removed bursts 
and selected for phases betweem -0.16 and -1-0.09. Firstly, 
we fitted the XIS and PIN dip spectra with a model of 
the form AG*(pcf*CUT -|- AGi*bb) i.e. applying the cov- 
ering factor PCF to the cut-off power law modelling the 
extended emission of the ADC, while the point-like black- 
body is subject to a column density AGi in addition to the 
Galactic column density AG, frozen at the non-dip value. 
An assumption of this model is that emission parameters 
cannot vary in dipping, only absorption parameters and 
so all of the continuum emission parameters of the black- 
body and cut-off power law were frozen at their non-dip 
values. 

Acceptable fits were obtained, with x^/d.o.f. ~ 1 but 

with clear residuals indicating the presence of absorption 
features. The emission line seen at ~6.4 keV in non-dip 
emission could not be seen in the residuals, and attempts 
to fit such a line resulted in a value of zero for the line nor- 
malization. Consequently no emission line was included in 
dip fitting. Inspection of the residuals revealed absorption 
lines at ^ 6.5, ^ 6.7, ^ 6.9 and '--^7.6 kcV. Corresponding 
lines were added to the spectral model. The line at 7.6 
keV appeared broader than the other lines in the residu- 
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Table 2. Best simultaneous fit to the non-dip and dip XIS -|- PIN spectra. 



Model Parameters Non-dip Dip 



JVh (10 atom cm ^) 






3.2±0.1 


3.2±0.1 


CUTOFF POWER LAW 
AiVn (10^^ atom cm"^) 
Covering fraction / 
photon index F 
cut-off energy (keV) 
normalization (lO"'^ ph cm" 


-2 s- 




0.0 
0.0 

J--O'_0.03 

nf--t-77 

'^■^-0.3 


22.2±0.8 
0.63±0.02 


BLACKBODY 

AATh (10^^ atom cm~^) 
temperature kTsB (keV) 
normalization (10^^ erg s~^ 


at 10 kpc) 


0.0 

1.35 ± 0.36 
8.3±5.1 


> 400 


ABSORPTION LINE 

energy E (keV) 
width a (kcV) 
normalization (lO"'' ph cm" 
equivalent width EW (eV) 


^ s 


) 




6.52±0.10 
0.01 (frozen) 
-1.0±0.8 
-18±15 


ABSORPTION LINE 
energy E (kcV) 
width a (keV) 
normalization (10~® ph cm" 
equivalent width EW (eV) 


-2 s- 




6.67±0.05 
0.01 (frozen) 
-1.24±0.71 
-12±7 


6.68tr,.;^ 

0.01 (frozen) 

-U.»_Q 9 

-15±15 


ABSORPTION LINE 

energy E (keV) 
width a (keV) 

normalization (10~' ph cm" 
equivalent width EW (eV) 


' s- 






6.94±0.14 
0.01 (frozen) 
-1.5±1.4 
-26±25 


ABSORPTION LINE 

energy E (keV) 
width (7 (keV) 
normalization (10^^ ph cm" 
equivalent width EW (eV) 


-2 s- 






7.6±0.2 
0.01 (frozen) 
-2.8±2.1 
-68±52 


dip/non-dip intensity shift 








0.79±0.04 


PIN normalization factor: 






1.12 ± 0.07 


1.12 (frozen) 


/d.o.f. 

X? 






1381/1634 
0.85 


140/263 
0.53 



90% confidence errors are given for all quantities except for the weaker 6.52 
and 6.68 keV lines where 67% confidence errors are shown for all parameters. 



als with a full width at zero height of ~0.6 kcV and so its 
width a was initially set at a corresponding 0.15 keV in 
the fitting, while the widths a of the first three lines were 
frozen at a value appropriate to narrow lines (0.01 keV). 

A formally acceptable fit could be obtained in this way. 

However, a comparison of the non-dip and deep dip spec- 
tra shown in Fig. 6 shows a vertically downwards shift in 
the dip spectrum, i.e. an energy-independent decrease of 



intensity. We investigated this, and determined the level of 
any energy-independent change during dipping by adding 
a multiplying factor to the spectral model, this factor af- 
fecting both cut-off power law and blackbody terms. All 
emission parameters were fixed at the non-dip values and 
the XIS-PIN normalization factor was fixed at the non-dip 
value of 1.12. There was a substantial improvement in the 
quality of fit with x^/d-O-f. becoming 140/263 compared 
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Suzaku MS + PIN XB 1323-619 




1 10 
Energy (keV) 

Fig. 6. Comparison of non-dip and deep dip spectra from best 
fits to XIS and PIN simultaneously in each. For clarity indi- 
vidual model components are not shown to demonstrate the 
energy-dependent nature of the dipping 



with the value without the vertical shift factor (164/263). 
Best-fit values for the final fit to the dip spectrum are 
shown in Table 2. The dipping was well explained by an in- 
crease in column density (ANh of at least 400x10'^^ atom 
cm~^ and covering of 63% of the Comptonized emission 
by absorber, with a smaller increase of TVh of 22 x 10^^ 
atom cm~^. This suggests that the blackbody is covered 
by central, denser regions of the absorbing bulge whereas 
the Comptonized emission is overlapped by lower density 
regions on average. Note that the best fit requires an un- 
absorbed fraction of 37%, i.e. part of the spectrum not 
overlapped is unabsorbed. In addition to the neutral ab- 
sorption, the dip spectrum requires an energy-independent 
decrease by a factor of 0.79±0.04. It was found that the 
covering factor / for the dip spectrum of 0.63±0.02 agreed 
well with the depth of dipping in the lightcurve of ~65%. 
Fig. 6 shows that below ~20 keV dipping is strongly 
energy-dependent as modelled by neutral absorber in our 
fitting. Above 20 keV the small, approximately energy- 
independent vertical shift can be seen. 

Table 2 also shows the best-fit line energies, all ob- 
tained by free fitting. The lines at 6.94 keV and 7.6 keV 
were detected at the 3-a level or higher, however, the other 
two lines were weaker and were only detected at l-cr and 
so we show these lines in the Table with l-cr errors. Fig. 
5 (lower panel) shows the absorption features using the 
best-fit continuum model with the lines removed from the 
spectral model. The line at ~6.5 keV apparently corre- 
sponds to the weak feature seen in the non-dip spectrum 
which we suggested was due to a less ionized Fc state. 
Possibly more than a single line is present, i.e. a forest of 
lines from Fe XVII onwards. There appears no doubt 



that the line at 6.68 keV is Fe XXV as seen in non-dip. The 

equivalent widths of these linos are greater than in non- 
dip; however, this may be due to the continuum decrease 
(EW = line flux/continuum intensity). But because of the 
large uncertainties in EW values, we can only say that 
the increase in EW is consistent with the change in con- 
tinuum. Indeed, it would not be expected that the bulge 
in the outer disk would cause increased absorption in high 
ionization states. The line at 6.94 keV, also seen in the 
observation with XMM (Church et al. 2005), can be iden- 
tified as Fe XXVI K^, and the line at 7.6 keV can be Ni 
XXVI. 

As a check, we also examined the energy-independent 
spectral change in the PIN spectrum alone, so as to be 
independent of XIS-PIN normalization using the energy 
band 20 - 70 keV where the effects of photoelectric ab- 
sorption are small. We carried out spectral fitting in this 
band using a cut-off power law (as the blackbody flux is 
negligible) fitting the non-dip spectrum, and then fixing 
the power law index and cut-off energy in fitting the dip 
spectrum. The normalization decreased by 24±3 %, agree- 
ing with the value of 21 ±4 % obtained from simultaneous 
fitting of XIS and PIN; however this result may be bi- 
assed by points at 20 keV where absorption makes a small 
contribution and we conclude that our best value for the 
degree by which electron scattering reduces the intensity is 
the factor shown in Table 3 from fitting both instruments 
simultaneously. 

4. Discussion 

Our spectral fitting results show that dipping consists 
of a large increase in photoelectric absorption for the 
Comptonized emission of the ADC, the intrinsic column 
density increasing from zero to 22 x 10^^ atom cm~^ in 
deep dipping, this being consistent with covering of an 
extended ADC source by an extended absorbing bulge in 
the outer disk as previously found for the dipping sources 
in general, and for XB 1323-619 in the observations with 
BeppoSAX (Baiucinska-Church et al. 1999), with RXTE 
(Barnard et al. (2001) and XMM (Church et al. 2005). The 
neutron star blackbody is subject to a high column den- 
sity which quickly removes this component in dipping as 
expected for a point-source and consistent with the source 
being covered by the highest column density part of the 
bulge in the outer disk, rather than an integration across 
the whole bulge as for the extended Comptonized emis- 
sion. 

An energy- independent shift can be seen at energies 
above 20 keV where photoelectric absorption is negligible 
indicating a degree of electron scattering, the shift being 
21±4%. A decrease of 10±10% was found in our obser- 
vation of the source with BeppoSAX in the band 20 - 50 
keV (Balucihska-Church et al. 1999). We can show that 
the actual high energy shift is consistent with the degree 
of electron scattering expected for the electrons present 
in the bulge in the outer disk. At an energy of 5 keV, 
the photoelectric absorption cross section ctpe is about 10 
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times the Thomson cross section <tt and the column den- 
sity for the Comptonizcd emission in dipping of 22 x 10^^ 
atom cm~^ reduces the intensity by a factor exp — A^h fPE 
~0.67 for the 63% of the emission that is covered, so that 
the overall factor is 0.67x0.63 + 0.37, i.e. a decrease of 
about 50% at 5 keV close to the actual decrease as seen 
in Fig. 6. 

At 20 kcV, (TpE is 10 times less than ctt and at higher 
energies becomes negligible. The electron density in the 
absorbing bulge in the outer disk will be between 1.0 and 
1.2 times the ion density depending on the exact ionization 
state, because of multiple electrons released in ionization 
of the metals so that for an electron column density Ne of 
^22x10^^ atom cm~^ electron scattering reduces the X- 
ray intensity of the covered part of the extended emission 
by a factor 0.86, so that the overall factor is 0.86x0.63 -|- 

0. 37, i.e. 0.90, a 10% reduction in intensity, somewhat less 
than our measured decrease of 21%, although the errors 
are substantial. Our spectral fitting was carried out using 
the cross sections for neutral absorber. We show below 
that the ionization parameter ^ = L/nr^, where L is the 
luminosity, n is the plasma density, and r the radial po- 
sition, can be no more than 50 in the outer disk because 
of the large distance from the neutron star. For ^ smaller 
than ~15, the degree of ionization in the bulge will be low, 

1. e. there will be a large neutral atom component. For ^ 
between 15 - 50, only a low ionization state is achieved; for 
example, the oxygen ion may reach a state O II - O VI, but 
not the higher states O VIII or O IX (Kallman & McCray 
1982). The absorption cross sections are reduced from the 
neutral absorber values and in oxygen, ctpe is reduced by 
up to a factor of -2 for O V / O VI (Vcrncr et al. 1996). 
Thus spectral fitting would require a column density up 
to two times higher than for neutral absorber so that A^e 
would also be higher and the reduction in intensity by 
electron scattering would rise to 20%. At higher energies, 
e.g. 20 - 70 keV, the Klein-Nishina formula shows that 
the mean Thomson cross section in this band decreases by 
12%, but this has little effect on the reduction in intensity 
so that, the expected decrease would again be 20%. 

Thus application of the progressive covering absorp- 
tion model to the Extended ADC LMXB emission model 
provides a very good description of the spectral evolu- 
tion in dipping, including both the effects of photoelec- 
tric absorption and electron scattering in the bulge in the 
outer disk. Assuming that the bulge has a radial extent of 
about 10% of the disk radius 3.5 x 10^'^ cm) and using 
the measured column density in dipping for the extended 
emitter (Table 3) we can estimate an average density in 
the bulge of ^ 7 x lO^'' cm^"*. From this it follows that 
the ionization parameter ^ cannot be larger than ~50, so 
that the matter is in a relatively low ionization state 
~ 1000 is needed for a high ionization state). If our mea- 
sured column densities are too low because we employ a 
neutral absorber model, ^ would be even smaller. 

Suggestions have been made by Diaz-Trigo et al. (2006) 
and Boirin et al. (2005) that dipping in all of the dipping 
LMXB, i.e. the large decreases in continuum flux, can be 



explained not in terms of absorption in the bulge in the 

outer disk but in an ionized absorber with high ionization 
state (log^ > 3) suggesting that the absorber is located 
much closer to the neutron star, although they do not 
specify the radial position of the absorber. It is proposed 
from this fitting that a progressive covering description 
of dipping is not needed so that there is no need for an 
extended Comptonizing region. However, the questions of 
why the bulge in the outer disk does not cause strong 
absorption, and what azimuthal structure causes dipping 
in their model are not addressed. 

Detailed recent work does reveal the presence of highly 
ionized line features in the dipping LMXB. For example, 
Jimenez-Garatc et al. (2003) investigated such features 
in XBT 0748-676 and concluded that dipping was due to 
neutral abosrber in the outer disk plus some contribution 
from ionized absorber also located in the outer disk but 
in layers above the disk. Parmar et al. (2002) report de- 
tection of broad emission and narrow absorption features 
in X 1624-490 using XMM-Newton. Lines present in dip- 
ping indicated that the additional absorbing material had 
lower ionization state than that producing the narrow ab- 
sorption features seen in persistent emission. 

However, the idea of a lower ionization state in dip- 
ping was taken by Diaz-Trigo et al. (2005) to mean that a 
single region of high but changing ionization state could 
describe both non-dip and dip spectra, without refer- 
ence to the bulge in the outer disk. Spectral evolution 
in dipping in several sources was modelled using a model 
ABS*xabs*(bb-|-pl) where ABS is neutral absorber and 
XABS is ionized absorber. This model had free parameters 
^, ionized absorber column density iVion and a line width 
a. The main parameter changes in dipping were a decrease 
of log^ from 3.5 typically to 3.0, an increase of TVion by 
ten times, e.g. from 4 x 10^^ to 40 x 10^^ atom cm~^ and a 
small increase in the column density of neutral absorber, 
e.g. typically by up to a factor of two from 0.3 x 10^^ to 
0.7 X 10^^ atom cm~^. In essence, this fitting still explains 
dipping as photoelectric absorption but in a region where 
the net effective cross section of the mixture of elements is 
reduced as the metals are in higher ionization states. The 
continuum change in dipping can be fitted, but the high ^ 
has the consequence that the absorber has to be located 
much closer to the neutron star than the bulge. 

Thus Diaz-Trigo et al. (2006) based on this spectral 
fitting, and Boirin et al. (2005), conclude that dipping 
can be explained without progressive covering of an ex- 
tended emission region as applied extensively to the dip- 
ping sources, e.g. by Church et al. (1997). They further 
suggest that an extended ADC is not required in their 
fitting, which is the basis of the Extended ADC emission 
model for LMXB in general (Church & Balucinska-Church 
2004) for which substantial evidence exists. Firstly, fitting 
non-dip and dip spectra in any dipping LMXB shows that 
the majority of the X-ray flux, i.e. the non-thermal emis- 
sion, is removed gradually and systematically proving that 
the emitter is extended (e.g. Church et al. 1997). The size 
of the Comptonization emission site has been measured by 
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dip ingress timing for all of the dipping LMXB (Church & 

Balucinska-Church 2004) showing that this emitter is in- 
deed very large, having a radial extent typically 50 000 km 
or 5 - 50% of the accretion disk radius. Strong indepen- 
dent evidence for an extended hot ADC is given by Schulz 
et al. (2009) who detect a variety of broad emission lines 
of highly ionized species including NeX La, MgXII L^, 
Si XIV L„, SXVI La, FcXXIV, FcXXV and FcXXVI in 
Chandra grating spectra of CygnusX-2. The emissivities 
indicate log^ > 3, and the radial positions of the emitters 
were found to be between 20 000 and 110000 km from the 
neutron star. It was thus concluded that these lines orig- 
inate in a hot, extended ADC and the radial distances 
arc in excellent agreement with the results of dip ingress 
timing. Schulz et al. (2008) in a study of CirX-1 similarly 
found a large number of emission line features rich in H- 
like and He-like features of high Z-clcmcnts such as Si, S, 
Ar and Ca. The spectra could be modelled as emission 
from a hot (T ~ 10^ K) photoionized plasma with ^ = 
1000, located at a radial distance r 10 000 km from the 
neutron star, i.e. a hot extended accretion disk corona. 

Given the now apparent overwhelming evidence for the 
existence of extended ADCs, it now seems necessary to in- 
clude such an extended emitter in any model fitting (i.e. 
by a progressive covering term), which is quite in contrast 
to simply using an ionized absorber as done, for example, 
in Diaz-Trigo et al. (2006). In addition, there are a number 
of elements lacking in the model which appear difficult to 
explain in terms of the model, as follows, which arc not 
a problem for the standard explanation of dipping. It is 
not specified at what radial position the ionized absorber 
exists, nor is it suggested what caiiscs the azimiithal struc- 
ture required in the absorber to cause dipping in a partic- 
ular range of orbital phases. It is difficult to see why there 
should be any azimuthal structure in an inner region where 
^ is high, especially since such structure would have to be 
similar to the known structure in the outer disk, i.e. the 
bulge. Finally the model docs not explain why the bulge 
should have no absorbing effect contrary to accepted ideas. 
Contrasted with this should be the existing explanations 
of dipping as progressive covering of the neutron star and 
an extended ADC by the extended absorbing bulge, which 
is physically reasonable, with substantial evidence for it. 

From the continuum changes, we turn to the line fea- 
tures detected in the Suzaku observation (Table 2). A pos- 
sible emission line consistent with 6.4 keV was seen in 
the Suzaku non-dip spectrum which was absent in the dip 
spectra. A similar feature had been seen in the 1997 RXTE 
spectrum at 6.43±0.21 keV (Barnard et al. 2001) although 
the line was not found in the 2003 XMM observation when 
the source was much brighter. Absorption lines were also 
seen: in the non-dip spectrum the only line we can be 
confident of was the 6.67 keV Fe XXV line. In dipping, 
four absorption lines were seen as discussed previously. 
The general presence of these medium to high ionization 
features are consistent with formation in a hot ADC and 
this hot extended ADC is, of course, a major element of 
the "Extended ADC" continuum emission model (Church 



& Balucinska-Church 2004). Other authors, in particular 

Schulz ct al. (2008, 2009) have presented strong evidence 
that similar lines originate in an extended, hot ADC. 

In the Suzaku observation of XB 1323-619 the luminos- 
ity was at a low level in contrast with the highest luminos- 
ity exhibited by the source in the 2003 XMM observation, 
after which the luminosity but reduced by about a factor 
of two in the RXTE observation of 2003, and fell further 
by almost another factor of two in the present observation. 
The X-ray burst rate has also decreased substantially to 
about the value it had in the Exosat observation, as shown 
in Fig. 2. The present burst rate supports the idea that 
the burst rate follows a definite burst rate-luminosity re- 
lation which has not changed over 20 years. We found in 
the present observation that the typical long burst de- 
cays having two time constants indicative of burning by 
the rp-proccss. Such behaviour was found previously in 
the source prior to the XMM observation. However, in the 
XMM observation it wasd found that the long tails were 
not seen in bursts indicating less rp-burning (Balucinska- 
Church et al. 2008), a possible reason being come complete 
burning of H to He between bursts at the higher luminosity 
and mass accretion rate. The return of the burst properties 
of the source to those determined for the Exosat, RXTE 
and BeppoSAX observations prior to the XMM observa- 
tion support the idea that the burst behaviour is deter- 
mined primarily by the mass accretion rate. 

The evidence is that the luminosity measured in the 
Exosat and Suzaku observations is a lower limit for the 
source. Between 1985 and 2003, the luminosity increased 
by a factor of about 4, implying that the mass accretion 
rate inc;reased by that factor. We cannot say whether this 
behaviour is periodic or not, however, the observed grad- 
ual increase over an 18 year period followed by a much 
faster decrease in luminosity over 4 years does not sug- 
gest a periodic change, and further observations will reveal 
whether the luminosity will increase as it did from 1985 
onwards. If the source were periodic, it could possibly be 
a triple system, in which case the third body would need 
to have a period of the order of 30 days. 
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